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After reviewing the modified Newtonian dynamics (MOND) proposal, we advocate that the 
associated phenomenology may actually not result from a modification of Newtonian gravity, 
but from a mechanism of "gravitational polarization" of some dipolar medium playing the 
O ' role of dark matter. We then build a relativistic model within standard general relativity to 

describe (at some phenomenological level) the dipolar dark matter polarizable in a gravita- 
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t/^ . tional field. The model naturally involves a cosmological constant, and is shown to reduce 

i ' to the concordance cosmological scenario (A-CDM) at early cosmological times. From the 

mechanism of gravitational polarization, we recover the phenomenology of MOND in a typi- 
cal galaxy at low redshift. Furthermore, we show that the cosmological constant A scales like 

^ ' flo, where ao is the constant MOND acceleration scale, in good agreement with observations. 

o 

^ I , 1 Introduction 

^— V . The mysteries of the nature of dark matter and dark energy are perhaps the most important 

OO ! ones of contemporary cosmology. Dark matter, which accounts for the observed discrepancy 

^^ ' between the dynamical and luminous masses of bounded astrophysical systems, is usually for- 

mulated within the so-called particle dark matter approach, in which the dark matter consists 
^ , of unknown non-baryonic particles, e.g. neutralinos as predicted by super-symmetric extensions 

$H ' of the standard model of particle physics (see^ for a review). Furthermore, the dark matter 

triggers the formation of large-scale structures by gravitational collapse and explains the distri- 
bution of baryonic matter from galaxy cluster scales up to cosmological scales by the non-linear 
growth of initial perturbations. Simulations suggest some universal dark matter density profile 
around distributions of ordinary baryonic matter l^. An important characteristic of dark mat- 
ter, required by the necessity of clustering matter on small scales, is that it should be cold (or 
non-relativistic) at the epoch of galaxy formation. Together with the observational evidence 
of dark energy (presumably a cosmological constant A) measured from the Hubble diagram of 
supernovas, the particle dark matter hypothesis has yielded the successful concordance model 
of cosmology called A-CDM, which reproduces extremely well the observed cosmic microwave 
background spectrum'^. 

However, despite these successes at cosmological scales, the particle dark matter approach 
has some difficulties'^ at explaining in a natural way the flat rotation curves of galaxies, one of 
the most persuasive evidence for the existence of dark matter, and the Tully-Fisher empirical 
relation between the observed luminosity and the asymptotic rotation velocity of spiral galaxies. 
In order to deal with these difficulties, all linked with the properties of dark matter at galactic 
scales, an alternative paradigm has emerged in the name of the modified Newtonian dynamics 



(MOND)I I I I Although MOND in its original formulation cannot be considered as a viable 
physical model, it is now generally admitted that it does capture in a very simple and powerful 
"phenomenological recipe" a large number of observational facts, that any pertinent model of 
dark matter should explain. 

It is frustrating that the two alternatives A-CDM and MOND, which are successful in com- 
plementary domains of validity (say the cosmological scale for A-CDM and the galactic scale 
for MOND), seem to be fundamentally incompatible. In the present paper, we shall propose a 
different approach, together with a new interpretation of the phenomenology of MOND, which 
has the potential of bringing together A-CDM and MOND into a single unifying relativistic 
model for dark matter and dark energy. This relativistic model will be shown to benefit from 
both the successes of A-CDM at cosmological scales, and MOND at galactic scales. 

2 The modified Newtonian dynamics (MOND) 

The original idea behind MOND'^'^^ is that there is no dark matter, and we witness a violation 
of the fundamental law of gravity (or of inertia). MOND is designed to account for the basic 
features of galactic dark matter halos. It states that the "true" gravitational field g experienced 
by ordinary matter, say a test particle whose acceleration would thus he a = g, is not the 
Newtonian gravitational field g^, but is actually related to it by 

^^{^)9 = m. (1) 

Here /i is a function of the dimensionless ratio g/ao between the norm of the gravitational field 
g = \g\, and the constant MOND acceleration scale ao — 1-2 x 10^^*^ m/s^, whose numerical 
value is chosen to fit the data. The specific MOND regime corresponds to the we ak gr avitv 
limit, much weaker than the scale oq. In this regime (where formally g ^ 0) we 
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J^)=^ + 0{g'). (2) 

\aoJ ao 

On the other hand, when g is much larger than oq (formally g -^ -l-oo) the usual Newtonian law 
is recovered, i.e. /i — > 1. Various functions /i interpolating between the MOND regime and the 
Newtonian limit are possible, but most of them appear to be rather ad hoc. Taken for granted, 
the MOND "recipe" ((l|)-(l2]) beautifully predicts a Tully-Fisher relation and is very successful at 
fitting the d etai led shape of rotation curves of galaxies from the observed distribution of stars 
and gas (seel^^for reviews). So MOND appears to be more than a simple recipe and may well 
be related to some new fundamental physics. In any case the agreement of ([I])-® with a large 
number of observations calls for a clear physical explanation. 

Taking the divergence of both sides of ([T]), and using the usual Poisson equation for the 
Newtonian field g^, we obtain a local formulation of MOND in the form of the modified Poisson 
equation^Liil 

V-(/iff) = -47rG/>b, (3) 

where pb is the density of baryonic matter. This equation can be derived from a Lagrangian, and 
that Lagrangian has been the starting point for constructing relativistic extensions of MOND. 
Such extensions postulate the existence of extra (supposedly fundamental) fields associated with 
gravity besides the spin-2 field of general relativity. Promoting the Newtonian potential U to 
a scalar field (j), scalar-tensor theories for MOND have been constructed ^^^ but shown to be 
non viable; essentially because light signals do not feel the presence of the scalar field, since 
the physical (Jordan-frame) metric is conformally related to the Einstein-frame metric, and the 



Maxwell equations are conformally invariant. This is contrary to observations: huge amounts of 
dark matter are indeed observed by gravitational (weak and strong) lensing. 

Relativistic extensions of MOND that pass the problem of light deflection by galaxy clusters 
have been shown to require the existence of a time- like vector field. The prototype of such the- 
ories is the tensor-vector-scalar (TeVeS) theory I '^ l I, whose non-relativistic limit reproduces 
MOND, and which has been extensively investigated in cosmology and at the intermediate scale 
of galaxy clusters ^^. Mo dified gra vitv theories such as TeVeS have evolved recently toward 
Einstein-aether like theories^ ^ ^ I Still, recovering the level of agreement of the A-CDM sce- 
nario with observations at cosmological scales remains an issue for such theories. In the present 
paper we shall follow a completely different route from that of modified gravity and/or Einstein- 
aether theories. We shall propose an alternative to these theories in the form of a specific modified 
matter theory based on an elementary interpretation of the MOND equation ([3]) . 

3 Interpretation of the phenomenology of MOND 

The physical motivation behind our approach is the striking (and presumably deep) analogy 
between MOND and the electrostatics of dielectric media^^. From electromagnetism in dielectric 
media we know that the Maxwell-Gauss equation can be written in the two equivalent forms 

V-E = — (cTfree + CTpol) ^=^ V • (Sr -E) = — CTfree , (4) 

where E is the electric field, afree and Upoi are the densities of free and polarized (electric) 
charges respectively, and Er = 1 + Xe is the relative permittivity of the dielectric medium. Such 
an equivalence is only possible because the density of polarized charges reads cjpoi = — V-P, 
where the polarization field P is aligned with the electric field according to -P = £oXe E, the 
proportionality coefficient Xc{E) being known as the electric susceptibility. 

By full analogy, one can write the MOND equation ([3]) in the form of the usual Poisson 
equation but sourced by some additional distribution of "polarized gravitational masses" ppoi 
(to be interpreted as dark matter, or a component of dark matter), namely 

V-g = -4^G(pb + /0poi) ^^ V-(/ig) = -4^Gpb- (5) 

This rewriting stands as long as the mass density of polarized masses appears as the divergence 
of a vector fied, namely takes the dipolar form 

Ppoi = -v-n, (6) 

where 11 denotes the (gravitational analogue of the) polarization field. It is aligned with the 
gravitational field g (i.e. the dipolar medium is polarized) according to 

Here the coefficient Xi which depends on the norm of the gravitational field g = \g\ in complete 
analogy with the electrostatics of dielectric media, is related to the MOND function by 

/i = l + X. (8) 

Obviously x can be interpreted as a "gravitational susceptibility" coefficient, while // itself can 
rightly be called a "digravitational" coefficient. It was shown in '—that in the gravitational case 
the sign of x should be negative, in perfect agreement with what MOND predicts; indeed, we 
have // < 1 in a straightforward interpolation between the MOND and Newtonian regimes, hence 



X < 0. This finding is in contrast to electromagnetism where Xe is positive. It can be viewed as 
some "anti-screening" of gravitational (baryonic) masses by polarization masses — the opposite 
effect of the usual screening of electric (free) charges by polarization charges. Such anti-screening 
mechanism results in an enhancement of the gravitational field a la MOND, and offers a very 
nice interpretation of the MOND phenomenology. Furthermore, it was pointed out that the 
stability of the dipolar dark matter medium requires the existence of some internal force, which 
turned out to be simply (in a crude quasi-Newtonian model ■''^) that of an harmonic oscillator. 
This force could then be interpreted as the restoring force in the gravitational analogue of a 
plasma oscillating at its natural plasma frequency. Finally, it seems from this discussion that the 
dielectric interpretation of MOND is deeper than a mere formal analogy. However the model^^ 
is clearly non-viable because it is non-relativistic, and it involves negative gravitational-type 
masses and therefore a violation of the equivalence principle at a fundamental level. 

4 Relativistic model for the dipolar dark fluid 

Here we shall take seriously the physical intuition that MOND has something to do with a 
mechanism of gravitational polarization. We shall build a fully relativistic model based on a 
matter action in standard general relativity. Note that this means we are changing the point 
of view of the original MOND proposal. Instead of requiring a modification of the laws of 
gravity in the absence of dark matter, we advocate that the phenomenology of MOND results 
from a physically well-motivated mechanism acting on a new type of dark matter, very exotic 
compared to standard particle dark matter. Thus, we are proposing a modification of the dark 
matter sector rather than a modification of gravity as in TeVeS like theories. 

4.1 Action and equations of motion 

From the previous discussion, the necessity of endowing dark matter with a new vector field 
to build the polarization field is clear. However this vector field will not be expected to be 
fundamental as in TeVeS like theories. Extending previous work ■ , our model will be based on 
a matter action (in Eulerian fluid formalism) in general relativity of the form 

S = JA''x./^gL[j^,eA^9,u]. (9) 

This action is to be added to the Einstein-Hilbert action for gravity, and to the actions of all 
the other matter fields. It contains two dynamical variables: (i) a conserved current J^ = au^ 
satisfying V^J'^ = 0, where m^ is the normalized four-velocity and a = {—J^J'^Y''^ is the 
rest mass energy density (we pose c = 1 throughout); (ii) the vector field ^^ representing the 
dipole four- vector moment carried by the fluid particles. This extra field being dynamical, the 
Lagrangian will also depend on its covariant derivative ^ui^] but in our model this dependence 
will occur only through the covariant time derivative £,^^ = u'^'V^S,^. The Lagrangian explicitly 
reads {see'^^ for details) 

1 



L = a 



-i-H + ^e.^- 



w(nx) , (10) 



where ^ = { (^/j — C^) (^^ ~ C^) } • The first term is a mass term in the ordinary sense (i.e. the 
Lagrangian of a pressureless perf ect fluid) , the second one is inspired by the action of spinning 

lool 

particles in general relativity i^^, and the third term is a kinetic term for the dipole moment. 
Finally, the last term represents a potential W describing some internal interaction, function 
of the polarization (scalar) field H^ = (J-^^jy H^H^)^'^, where H^ = a£^^ is the polarization 
four-vector, and -L^^, = g^y + n^Uj^ is the projector orthogonal to the four-velocity. 



By varying the action (|9j)- (jl0p with respect to both the current J^ and the dipole moment 
^^, we get two dynamical equations: an equation of motion for the dipolar fluid, and an evolution 
equation for the dipole moment ^^. From these equations it can be shownl^that we can impose 
the constraint H = 1 as a particular way of selecting a physically interesting solution, such that 
the final equations depend only on the space-like projection ^^ =_L^j^^'' of the dipole moment 
^^. The final equations we obtain are 



u^' = -T^ = -i^lW\ (11) 

1, 
a 



nf" = -v^ (w - n^W) - ClR^'p^x^^u^ > (12) 



where we denote O'^ = n^(l + ^±W) +±u ^'^, and employ the notations ^^ = S.'^/^x = n^/n_L 
and W' = dW/dllx. The motion of the dipolar fluid as given by (jlip is non-geodesic, and driven 
by the internal force J-'^ derived from the potential W. Observe the coupling to the Riemann 
curvature tensor in the equation of evolution ()12p of the dipole moment. By varying the action 
with respect to the metric gf^i, we obtain the stress-energy tensor T^'^. Using the canonical 
decomposition T^'^ = ru^u'^ + V 1-^'^ + 2Q^^u'^' + Ti^'^ , we find the energy density r, pressure 
V, heat flux Q^ (such that u^Q^ = 0) and anisotropic stresses Tj^'^ [uyT^^^^ = and S[^ = 0) as 

r = >V-n^>V'+p, (13a) 

p = ->V + ^nxW', (13b) 

gA* = CT ^^ + Ii^W'u^' - n^ Va-u^ , (13c) 

Here the contribution p to the energy density involves a monopolar term a and a dipolar term 
— VaII^ which clearly appears as a relativistic generalisation of ([6]), and will play the crucial 
role when recovering MOND: 

p = a- VaB^ . (14) 

4-2 Weak field expansion of the internal potential 

The dipolar fluid dynamics in a given background metric, and its influence on spacetime are now 
known; in the following we shall apply this model to large-scale cosmology and to galactic halos. 
For both applications we shall need to consider the model in a regime of weak gravity, which 
will be either flrst-order perturbations around a Friedman-Lemaitre-Robertson- Walker (FLRW) 
background in cosmology, or the non-relativistic limit for galaxies. A crucial assumption we 
make is that the potential function W admits a minimum when the polarization 11^ is zero, and 
can be Taylor-expanded around that minimum, with coefflcients being entirely specifled (modulo 
an overall factor G) by the single surface density scale built from the MOND acceleration ao. 

These coefficients in the expansion of W when n_L — > will be fine-tuned in order to recover the 
relevant physics at cosmological and galactic scales. Physically, this expansion corresponds to 
n^ <C S and is valid in the weak gravity limit g <C ao. Clearly, the minimum of W is nothing 
but a cosmological constant A, and we find 

>v(n^) = ^ + 27rG ni + ^ ni + o(nl) . (i6) 

The expansion is thereby determined up to third order inclusively. 



Our assumption that the function W involves the single fundamental scale S implies in 
particular that the cosmological constant A should scale with G'^T? ~ Oq. We thus introduce a 
dimensionless parameter a through 

aao = — A/- = aA. (17) 

This parameter represents a conversion factor between oq and the natural acceleration scale ca 
associated with the cosmological constant A. Posing x = n^/E, we find that (J16p can be recast 
in the form W(n_|_) = GvrGS^ w{x), where 

w{x) = o?T? + -x2 + -x3 + Oix^) . (18) 

o 9 

The present model should be considered only as "effective" or phenomenological, in the sense 
that the weak-gravity expansion (J16p or (jlSp should come from a more fundamental (presum- 
ably quantum) underlying theory. Therefore, the coefficients in (jlSp should not be given by 
exceedingly large or small numbers (like lO^'' or 10"^'^), but rather be numerically of the order 
of one, up to a factor of, say, ten. Hence, we expect that a itself should be around one, and 
we see by (J17p that the cosmological constant A should naturally be numerically of the order 
of Oq. Notice though that our model does not provide a way to compute the exact value of a. 
However, one can say that the "cosmic coincidence" (see e.g.l^ between the values of oq and 
OA — with the measured conversion factor being a ~ 0.8 — finds a natural explanation if dark 
matter is made of a fluid of polarizable gravitational dipole moments. 

5 Recovering the A-CDM scenario at cosmological scales 

Consider a small perturbation of a background FLRW metric valid between, say, the end of 
the inflationary era and the recombination. The dipolar fluid is described by its four- velocity 
u^ — ui^-\-^u^ ^ where bu^ is a perturbation of the background comoving four-velocity vP" = (1, 0), 
and by its rest mass density o" = ct -|- ba^ where ba is a perturbation of the mean cosmological 
value a. The crucial point in our analysis of the dipolar fluid in cosmology is that the background 
value of the dipole moment fleld .^^ (which is orthogonal to the four-velocity and therefore is 
space-like) must vanish in order to preserve the spatial isotropy of the FLRW background. We 
shall thus write ^j^ = b^j_, and similarly H^ = bH^ for the polarization. 

At flrst perturbation order, the stress-energy tensor with explicit components (J13p can be 
naturally recast, using also ([TH]) . in the form T^'^ = Tl^^ + ^dm; where the explicit expressions 
of the dark energy and dipolar dark matter components are 
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^d: = -^5'^^ (19) 



TZ = puf'u'' + 2QK''K (20) 

At that order, we flnd that the dipolar dark matter density reduces to p given by ()14p . The heat 
flux in ()20p reads as Q^ = <^ C± ~ H^Va^'^; it is non zero, however it is perturbative because 
so are both ,^j^ and H^^. The point for our purpose is that at linear perturbation order Q^ 
can be absorbed into a redefinition of the perturbed four-velocity of the dipolar fiuid. Posing 
bu^ = bu^ + Q^/W and introducing the effective four- velocity u^ = u^ + bu^, we find that at 
flrst order the dipolar dark matter fluid is described by the stress-energy tensor 

TZ = P^'^''^ (21) 



which is that of a perfect fluid with four- velocity u^, vanishing pressure and energy density ([li 
In the linear cosmological regime, the dipolar fluid therefore behaves as standard cold dark 



matter (a pressureless fluid) plus standard dark energy (a cosmological constant). Adjusting the 
background value a so that 0,dm — 23%, the model is thus consistent with the standard A-CDM 
scenario and the cosmological observations of the CMB fluctuations (seel^for more details). 

6 Recovering the MOND phenomenology at galactic scales 

Next, we turn to the study of the dipolar dark fluid in a typical galaxy at low redshift. We 
have to consider the non-relativistic limit (c -^ +00) of the model; we consistently neglect all 
relativistic terms 0(c~^). It is straightforward to check that in the non-relativistic limit the 
equation of motion (jlip of the dipolar fluid reduces to 

^=9- n^W , (22) 

where g is the local gravitational field generated by both the baryonic matter in the galaxy and 
the dipolar dark matter. Applying next the standard minimal coupling to gravity in general 
relativity, we find that the gravitational field obeys the Poisson equation 

V-g = -A7TG{py, + p), (23) 

where /Ob and p are respectively the baryonic and dipolar dark matter mass densities. From (|14p 
the dipolar dark matter density reduces in the non-relativistic limit to 

p = a- V-n^ . (24) 

The first term is a usual monopolar contribution: the rest mass density a of the fluid, while 
the second term is the dipolar contribution, and can be interpreted as coming from the fluid's 
internal energy. Here 11^ denotes the spatial components of the polarization. 

In order to recover MOND, we need two things: (i) to find a mechanism for the alignment 
of the polarization field 11^ with the gravitational field g, so that a relation similar to ([7]) 
will apply; (ii) to justify that the rest mass density a of dipole moments in (I24p is small with 
respect to the baryonic density pb, hence the galaxy will mostly appear as baryonic in MOND 
fits of rotation curves. We have proposed inl^^ a single mechanism able to answer positively 
these two points. We call it the hypothesis of weak clustering of dipolar dark matter; it is 
an hypothesis because it has been conjectured but not proved, and should be checked using 
numerical simulations. The weak clustering hypothesis is motivated by a solution of the full set 
of equations describing the non-relativistic motion of dipolar dark matter in a typical baryonic 
galaxy whose mass distribution /Ob is spherically symmetric (see the Appendix inl^^. This 
particular solution corresponds to an equilibrium configuration in spherical symmetry, for which 
V = and a = (To(r'). The dipole moments remain at rest because the gravitational field g is 
balanced by the internal force ^ = tl±W'; for that solution the right-hand-side of (j22p vanishes. 

During the cosmological evolution we expect that the dipolar medium will not cluster much 
because the internal force may balance part of the local gravitational field generated by an 
overdensity. The dipolar dark matter density contrast in a typical galaxy at low redshift should 
thus be small, at least smaller than in the standard A-CDM scenario. Hence o" <C /Ob; and we 
could even envisage that a stays around its mean cosmological value, o" ~ ct ^ pb. Now, because 
of its size and typical time-scale of evolution, a galaxy is almost unaffected by the cosmological 
expansion of the Universe. The cosmological mass density a of the dipolar dark matter is not 
only homogeneous, but also almost constant in this galaxy. The continuity equation reduces to 
V-(a'u) ~ 0, and the most simple solution corresponds to a static fluid verifying d ~ 0. By (j22p 
we see that the polarization field n_|_ is then aligned with the gravitational field g, namely 

g ~ U±W' . (25) 



On the other hand, because a ^ pb by the same mechanism, we observe that the gravitational 
field equation (I23p with (j24p becomes 

V-{g-A7:GU±)^-ATTGpb, (26) 

which according to ([5]) is found to be rigorously equivalent to the MOND equation. Finally, by 
inserting the expression of the potential (|16|) into (|25|) and comparing with the defining equation 
([7]), we readily find the MOND behaviour of the gravitational susceptibility coefficient as 

x(5) = -l + - + 0(<?'), (27) 

in complete agreement with ([2]). We can thus state that the dipolar fluid described by the action 
([9])- (jl0p explains the phenomenology of MOND in a typical galaxy. Note also that in this model 
there is no problem with the light deflection by galaxy clusters. Indeed the standard general 
relativistic coupling to gravity implies the usual formula for the bending of light. 

To conclude, the present model reconciles in some sense the observations of dark matter on 
cosmological scales, where the evidence is for cold dark matter, and on galactic scales, which 
is the realm of MOND. In addition, it offers a nice unification between the dark energy in the 
form of A and the dark matter in galactic halos. More work should be done to test the model, 
either by studying second-order perturbations in cosmology, or by computing numerically the 
non-linear growth of perturbations and comparing with large-scale structures, or by studying 
the intermediate scale of clusters of galaxies. 
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